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Based on theoretical analysis of PCM (Phase Change Material) solidification process, the model of improved void 


cavity distribution tending to high temperature region is established. Numerical results are compared with NASA 


(National Aeronautics and Space Administration) results. Analysis results show that the outer wall temperature, 


the melting ratio of PCM and the temperature gradient of PCM canister, have great difference in different void 


cavity distribution. The form of void distribution has a great effect on the process of phase change. Based on si- 


mulation results under the model of improved void cavity distribution, phase change heat transfer process in 


thermal storage container is analyzed. The main goal of the improved designing for PCM canister is to take 


measures in reducing the concentration distribution of void cavity by adding some foam metal into phase change 


material. 
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Introduction 


A heat pipe receiver with high temperature phase 
change material is a very important component for ad- 
vanced solar dynamic power generation system during 
traversal of low earth orbit spacecraft through the dark 
phase of their orbit cycles. The main functions of the 
receiver are the absorption of solar energy, thermal en- 
ergy storage (TES) for the dark period, and the heating of 
working gas [1-3]. 

A heat pipe receiver contains an array of heat pipes 
with mounted cylindrically toroidal thermal energy sto- 
rage modules. A PCM canister is a basic thermal storage 
unit of a heat pipe receiver. Void cavity appears in canis- 
ter when PCM freezes and shrinks. High temperature 
phase change material is commonly used as 80.5%LiF- 
19.5%CaF>. Its void ratio is 8%-22%. Since the thermal 
resistance of void cavity is very high, void cavity pre- 
vents phase change. Thermal spot and thermal ratcheting 
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appear near the concentrated void cavity easily. Stress 
and creep deformation will appear with the effect of the 
thermal-structural coupling of void cavity and high tem- 
perature heat load. Thus the usage life of high tempera- 
ture phase change thermal storage canister will decrease 
and even destroy. The thermodynamic performance of 
high temperature phase change material canister and 
compound phase change material canister need to be 
deeply studied. The influence of void cavity on thermal 
performance of PCM canister for heat pipe receiver un- 
der microgravity needs to be analyzed in detail. 


Structural Design and Working Principle of a 
Heat Pipe Receiver 


The heat pipe receiver design is shown in Fig. 1 [4]. 

The receiver is constructed in three axial sections, a 
receiver portion, a TSD (Thermal Storage Device) por- 
tion, and a HSHX (Heat Source Heat Exchanger) por- 
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Nomenclature 
A area (mô 
Almich the mushy zone constant 
C specific heat capacity at constant pressure 
P (J/ (kg K) 
H enthalpy (J/kg) 
AH latent heat of phase change (J/kg) 
h sensible enthalpy (J/kg) 
k thermal conductivity (W/(m: K)) 
L the axial length (m) 
Q heat flux (W) 
R radial length (m) 
Rco radial length of canister outer wall (m) 
Rinick radial thickness of PCM (m) 
r radial direction 


TSD Portion 


Working Fluid Inlet _-] 
“SHSHX Portion 


Heat Pipe Plate Fin Heat Exchanger 


Fig. 1 The configuration of the heat pipe receiver 


tion.The first of these sections is located nearest the re- 
ceiver aperture and includes no PCM canisters. Solar 
energy entering through the aperture strikes the central 
tubes directly. Inside the tube is a heat pipe working fluid. 
The evaporated fluid then transports the energy into the 
remaining two sections. In the central section the tube is 
surrounded by PCM canisters, with a configuration much 
like those in the base case design. The central section 
contains 91.6 kg of LiF»CaF divided into full length 
canisters (4.8-cm ID, 8.2-cm OD, .71-m length) sur- 
rounding each of the 20 tubes. 


Numerical Simulation of PCM Canister for Heat 
Pipe Receiver under Microgravity 


Physical Model 


The heat transfer of PCM canister of thermal storage 
unit of the heat pipe receiver [5] includes 

(1) Heat conduction among solid, liquid in PCM zone 
and saturated vapor in void cavity, heat conduction 
among PCM canister wall, fin and other enhanced heat 
additions, 
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S source term 

T temperature (K) 
time(s) 

Z axial direction 


Greek letters 
B volume expansion coefficient of liquid 
(1/K) 


“ a small number (0.001) to prevent division 


by zero 
ġ turbulence quantity 
p density (kg/m*) 
Subscripts 
hw heat-pipe wall 
Ref reference 


(2) Radiation among the interfaces of void cavity with 
different temperature, radiation among different PCM 
crystals, and 

(3) Evaporation and condensation heat transfer be- 
tween high- and low- temperature interface in void cavity. 

The above process of heat transfer is transient. If all 
the above are considered, it is very difficult to solve the 
problem. To find easy solutions and simplify the calcula- 
tions, some assumptions are made as follows when the 
heat transfer model of the PCM canister is built. 

(1)The temperature along heat pipe wall is uniform. 

(2)The initial void volume of thermal storage unit is 
fixed, and void cavity distributes in the region tending to 
high temperature during the process of solidification. 
Prior to the domestic and foreign simulation of the phase 
change process in PCM canister, the existence of void 
cavity is usually ignored, or void cavity is simply as- 
sumed to distribute in annual region near the outer wall 
of PCM canister (hereinafter referred as the model of 
simple void cavity distribution). In this paper, the distri- 
bution model of void cavity tending to high temperature 
region during the process of solidification is established 
by assuming that the volume of void cavity occupies 10% 
of the volume of PCM canister. The detailed method is as 
follows. The solidification process of liquid phase PCM 
in thermal storage container is firstly simulated in the 
condition of ignoring void cavity, and then, simulation 
results of temperature field and PCM solidification can 
be obtained. In the end, the model of void cavity distri- 
bution is established by selecting 10% of the volume of 
final solidification region as void cavity (hereinafter re- 
ferred as the model of improved void cavity distribution). 
The model for improved void cavity distribution is more 
reasonable than that of the previously neglecting void 
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cavity or simple void cavity distribution, and the calcula- 
tion process is relatively simpler. 

(3)Heat convection of liquid PCM is ignored. Natural 
convection disappears under microgravity. 

(4)The contact resistance between the PCM canister 
and heat pipe wall is ignored. 

(5)The surfaces of all void cavities are the ones of 
gray bodies with diffusion reflection. 


Mathematical Model 


Mathematical Model of PCM Zone [6] 


Energy Equation 
For solidification-melting problems, the energy equa- 
tion is written as 


0(pH = 
PR Yov H)=V(KVT)+S (1) 
The enthalpy of the PCM is defined as 
H =h+AH (2) 
The sensible enthalpy of the PCM is defined as 
T 
h=Hyep + | CpdT (3) 
Tref 
The liquid fraction of the PCM is defined as 
B=0 if PST. 
= Fes if Typ <T<Tys (4) 
Tigs Tas 
ß=1 if T >To 


The solution for temperature is essentially an iteration 
between the energy equation (Equation 1) and the liquid 
fraction equation (Equation 4). Directly using Equa- 
tion (4) to update the liquid fraction usually results in 
poor convergence of the energy equation. In FLUENT 
6.2, the method suggested by Voller and Swaminathan [7] 
is used to update the liquid fraction. 


Momentum Equation 


The enthalpy-porosity technique treats the mushy re- 
gion (partially solidified region) as a porous medium. 
The porosity in each cell is set equal to the liquid fraction 
in that cell. In fully solidified regions, the porosity is 
equal to zero, which extinguishes the velocities in these 
regions. The momentum sink due to the reduced porosity 
in the mushy zone takes the following form 

1- BY > 3 
S =P g Ana V-V) (5) 

The mushy zone constant measures the amplitude of 
the damping. The higher this value, the steeper the transi- 
tion of the velocity of the material to zero as it solidifies. 
Very large values may cause the solution to oscillate. 
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Turbulence Equation 

Sinks are added to all of the turbulence equations in 
the mushy and solidified zones to account for the pres- 
ence of solid matter. The sink term is very similar to the 
momentum sink term (Equation 5). 


S= = Anush? (6) 
é 


Mathematical model for zone of void cavity [8] 

Since the axial temperature gradient is very small, the 
temperature of void cavity is defined by radial steady 
heat transfer equation, 


22 (rt <0 (7) 


The thermal conductivity of vapor in void cavity can 
be calculated according to the kinetic theory of gas. 


k, =1.457 x10” JT, (8) 


An Approach to Solution [9] 


Computational fluid dynamics (CFD) software 
(FLUENT 6.2) is used to simulate the heat transfer of 
PCM canister of the heat pipe receiver. Software 
(GAMBIT 2.1) is used as a pre-processor. In FLUENT 
6.2, the separate solver is used to solve the equations. 
Finite control volume (FCV) method is used to solve the 
above equations. Thus, non-linear differential equations 
are linearized. Alternate direction iteration (ADI) method 
is used as a method of solution. In FLUENT 6.2, a solidi- 
fication-melting model upon the enthalpy-porosity me- 
thod is specially provided to deal with phase changes. 


Mesh Generation 

Fig. 2 shows the geometric model and mesh divisions 
of thermal storage unit for heat pipe receiver with three 
kinds of void cavity distribution relevantly. The regions 
of PCM zone, void cavity and canister wall are plotted 
with the triangle meshes. The spacing of mesh is 0.5 mm. 


Boundary and Initial Conditions 


Boundary Conditions 


z=0,AH =0 (9) 
z=L,AH =0 (10) 
r= Reo, AH =0 (11) 
oT, 
r = Rhy, Kh A > +Q(Thy)/A (12) 


The boundary condition of heat pipe wall (canister in- 
ner wall) is taken as periodic various heat flux. The ca- 
nister outer wall and side wall are adiabatic with outside. 
In the void cavity, thermal radiation is considered. 
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Outer wall 


Side wall 
Side wall 
PCM 
Inner wall including heat pipe wall 
(a) Ignoring void cavity 
_ Outer wall 
Side wall 
Inner wall including heat pipe wall 
(b) Simple void cavity distribution 
Outer wall Void cavity 
PCM 
Side wal 
Side wal 
Inner wall including heat pipe wall 
(c) Improved void cavity distribution 
Fig. 2 Geometric model and mesh divisions of thermal 


storage unit for heat pipe receiver with three kinds of 
void cavity distribution 


The driving heat flux from the heat pipe wall to PCM 
canister is kept constant when the PCM thickness of the 
canister varies. The value of heat flux is positive during 
sunlight period, and it is negative during eclipse period. 
The heat flux is the function of heat pipe wall tempera- 
ture during both sunlight and eclipse period. 


Initial Conditions 


The initial temperature is set to 950 K. PCM is solid. 
Each orbit cycle consists of 94.61 minutes. The physical 


properties of the canister wall and phase change material 
are constant. The eutectic salt 80.5%LiF-19.5%CaF> se- 
lected as the PCM in our simulation, with a melting point 
of 1040 K, and the cobalt-base super alloy Haynes 188 as 
the canister material. 

The temperature of the heat sink is 200 K. The height 
of orbit is 500 km. The pressure of low earth orbit is 10° 
-10” Torr. 


1020 


Temperature/K 
S 
© 
© 


heat pipe wall 
-— — PCM canister outer wall 


0 100 200 300 400 500 600 700 800 900 
Time Step/min 


(a) Ignoring void cavity 


APOE 


— heat pipe wall 
=-=- PCM canister outer wall 


0 
0 100 200 300 400 500 600 700 800 900 1000 
Time Step/min 


(b) Simple void cavity distribution 


980 4 —— heat pipe wall 
=-=- PCM canister outer wall 


0 200 400 600 800 1000 1200 1400 
Time Step/min 
(c) Improved void cavity distribution 


Fig. 3 Temperature fluctuations on heat-pipe wall and PCM- 
canister outer wall during orbital periods with three 
kinds of void cavity distribution 


Results and Discussion 


Comparative Analysis on Calculation Results under 
the Models of Different Void Cavity Distribution 


Fig. 3a, Fig. 3b and Fig. 3c show the temperature fluc- 
tuations on heat-pipe wall and PCM-canister outer wall 
during orbital periods by our simulation corresponding to 
the models of neglecting void cavity, simple void cavity 
distribution and improved void cavity distribution re- 
spectively. Fig. 4 and Fig. 5 show the fluctuations of 
temperature and liquid fraction in one steady period re- 
spectively. It can be seen that with the periodic variation 
of incident solar heat flux, heat pipe wall temperature 
attains the maximum at the end of sunlight, and attains 
the minimum at the end of eclipse. 
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Fig. 4 Temperature fluctuations of heat pipe wall and PCM 
canister outer wall in one steady period with three 
kinds of void cavity distribution 
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Fig.5 Liquid fraction in one steady period with three kinds of 
void cavity distribution 


By comparison of the temperature change curve and 
liquid phase fraction distribution of PCM canister corre- 
sponding to three kinds of void cavity distribution, it can 
be seen that the outer wall temperature and the melting 
ratio of PCM have great difference under the models of 
neglecting void cavity, simple void cavity distribution 
and improved void cavity distribution. The form of void 
cavity distribution has a great influence on the process of 
phase change and thermal storage. The melting rate and 
effective utilization of PCM for three void cavity distri- 
bution models are shown in Table 1. 


Table 1 The melting ratio and effective utilization ratio of 
PCM for three kinds of void cavity models 


Melting ratio of Effective utilization 


Model of void cavity PCM tio of PCM 
ratio o 


Ignoring void cavity 0.0612-0.97393 0.9127 
Simple void cavity 

0.03405-0.8015 0.76745 
distribution 
Improved void cavity 

0.15998-0.9609 0.80092 


distribution 


It can be seen that the effective utilization ratio of 
PCM is the minimum, the middle, and the maximum 
respectively with the model of neglecting void cavity, 
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simple void cavity distribution and improved void cavity 
distribution. 

By the above comparative analysis, it can be found 
that different void cavity distribution model has a great 
influence on simulation results for the process of phase 
change thermal storage. Therefore, the establishment of a 
reasonable void cavity model is a precondition to im- 
prove the accuracy of simulation results. In the process of 
simulation, a simple void cavity model is used to assume 
that void cavity distributes in a uniform annular region 
close to the outer wall of thermal storage container, 
which simplifies the determination of void cavity region 
and mesh generation of geometric model, but this kind of 
void cavity distribution model is too simple, and will 
affect the accuracy of simulation results. The assumption 
of ignoring void cavity simplifies the calculation process 
greatly, but which is obviously inconsistent with the ac- 
tual situation and makes the simulation results have great 
difference with the actual process of phase change ther- 
mal storage. The improved void cavity model is estab- 
lished on the basis of the theoretical analysis. The im- 
proved model tending to high temperature region is based 
on the actual solidification process. The model of im- 
proved void cavity distribution is more reasonable and 
reliable than that of ignoring void cavity and simple void 
cavity distribution. According to the main characteristics 
of phase change thermal storage during the above process, 
the improvement for the designing of PCM canister is 
necessary. The concentrated distribution of void cavity 
not only reduces thermal performance of PCM canister 
due to local heat resistance, but also generates thermal 
stress and influences structural reliability of thermal 
storage container due to large temperature gradient. 
Therefore, the main goal of the improved designing for 
thermal storage container is to take measures in reducing 
the concentration distribution of void cavity by adding 
some foam metal into phase change material. In addition, 
heat transfer in the container can be effectively enhanced 
by side wall, which needs to be fully considered in the 
structural designing of thermal storage container. 


Comparative analysis on our numerical results and 
NASA numerical ones 


In NASA project [10], the model of void cavity dis- 
tribution is built by the assumption as follows. Void cav- 
ity distributes along the outer wall of thermal storage 
canister. During eclipse periods, PCM in the TSD region 
freezes. PCM releases latent heat to heat pipe. At the 
same time, PCM shrinks and void cavity appears. Since 
the temperature of PCM close to canister inner wall de- 
creases more quickly than that of PCM close to canister 
outer wall, PCM close to canister inner wall freezes ear- 
lier. PCM in canister freezes gradually along the radial 
from inner wall to outer wall. Most of void cavities ap- 
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pear near the outer wall in canister. 

In order to make a better comparison, the model of 
void cavity distribution for our simulation is the same 
with that for NASA project. Fig. 6 shows a comparison 
between our numerical results and NASA numerical ones 
[10]. In the comparison, it can be seen that the corres- 
ponding temperature differences between our simulation 
and NASA simulation on both heat-pipe wall and PCM- 
canister outer wall are small. The temperature-change 
trend of our numerical results agrees with that of NASA 
numerical ones. The difference is that phase change tem- 
perature of PCM (LiF) for our simulation is set to 1117- 
1125 K, which is more in line with the actual melting 
process of LiF. In NASA model, the phase transition 
temperature of 1121 K is adopted. The temperature fluc- 
tuation both on heat-pipe wall and in PCM canister re- 
mains less than 16 K throughout a sunlight and dark cycle. 
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Fig. 6 Liquid fraction in one steady period with three kinds 
of void cavity distribution 
Conclusions 


With both void cavity and phase change considered, 
PCM canister of heat pipe receiver is numerically simu- 
lated under microgravity. Some conclusions are as fol- 
lows. 

(1) The physical model and mathematical model of 
phase change during the process of heat transfer in ther- 
mal storage container are established. Based on theoreti- 
cal analysis and calculation results of PCM canister dur- 
ing the process of solidification, the model of improved 
void distribution is established. 

(2) The phase change thermal storage process is nu- 
merically simulated under the model of ignoring void 
cavity distribution, simple void cavity distribution and 
improved void cavity distribution. The calculation results 
for three kinds of void distribution are compared and 
analyzed. The results show that the form of void cavity 
distribution has a great effect on calculation results. 

(3) Our simulation results are compared with NASA 
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ones under the model of simple void cavity distribution. 
The accuracy of calculation model and method is verified. 

(4) Based on simulation results under the model of 
improved void cavity distribution, phase change heat 
transfer process in thermal storage container is analyzed. 
The main goal of the improved designing for PCM can- 
ister is to take measures in reducing the concentration 
distribution of void cavity by adding some foam metal 
into phase change material. 
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